Phenomenological theory of magnetic anisotropy 
in thin films and multilayers 



A.N. Bogdanov 1 * , U.K. R6filer 2 # 
1 Max-Planck-Institut fiir Physik komplexer Systeme 
Nothnitzer Strafie 38, D-01187 Dresden, Germany 
2 Institut fiir Festkorper- und Werkstofforschung Dresden 
Postfach 270016, D-01171 Dresden, Germany 
(February 1, 2008) 

Within a phenomenological approach the density of induced uniaxial anisotropy in nanostruc- 
tures K(r) is treated as a physical field additional to the magnetization M(r). The equilibrium 
distributions of K(r) are formed under the influence of the surfaces and internal interactions. The 
functions K (r) and M(r) are calculated for a magnetic layer between two nonmagnetic spacers. It 
is shown that the transition from the phase with magnetization in plane to the perpendicular phase 
occurs continuously via an intermediate phase which is inhomogeneous across the layer-thickness. 
The theory explains experimentally observed thickness dependences of the effective anisotropy and 
gives a method to obtain values of the characteristic parameters for the theory from experimental 
data. It is shown that the separation of the induced anisotropy into volume and surface contribution 
is valid only for thick films. 



Magnetic anisotropy is one of the key properties exploited in ferromagnetic nanostructures. Reduced dimensionality 
may greatly alter the effective magnetic anisotropy of such structures. For layered systems, numerous investigations 
demonstrate that reduced thickness and complex interactions on surfaces and interfaces induce uniaxial magnetic 
anisotropy which, by orders of magnitude, exceeds the values of the intrinsic anisotropy in the corresponding bulk 
materials |l| . During the last decade a large variety of magnetic thin films and multilayers have been synthesized and 
investigated experimentally. They have already found applications in modern magnetoelectronics such as spin-valves 
or magnetoresistive heads Q . 

Recent first principle numerical calculations have proved the reliability of the quantum-mechanical microscopic 
theory in its applications to magnetic nanostructures m . They are offering considerable insight into the phenomena 
of induced anisotropy. However, these ab initio calculations are still unable to give a complete description of the 
magnetization structures and processes in real layered systems. Hence, our understanding of and control over the 
magnetic anisotropy of nanostructures is rather incomplete. Up to now, the analysis of experimental data on anisotropy 
effects in magnetic layered systems is mostly based on the effective volume anisotropy method introduced by Neel 
about fifty years ago . Within this approach the value of average induced uniaxial anisotropy in a magnetic layer 
of thickness d is given by an empirical ansatz 

K cS = K v + -± (1) 

where Ky is the volume anisotropy (per unit volume), and K$ is the surface (interface) contribution (per unit 
area). For many layered systems with reduced thickness, the function K(l/d) becomes strongly nonlinear which 
is at variance with eq. ([T|) | . Moreover, the general validity of the Neel theory has been questioned in 
a number of publications [fOlW. The separation of the uniaxial anisotropy into volume and surface contribution 
and the reduction of surface anisotropy into an effective volume contribution should be considered as an excessive 
simplification. Furthermore, equation ([!]) implies that the anisotropy is constant within a magnetic layer and, thus, 
stabilizes a homogeneous distribution of the magnetization. It was shown, however, that the competition between 
volume and surface anisotropies may induce inhomogeneous states in magnetic films and multilayer structures [fl3|| . 
The surface induced anisotropy can also gradually relax into the depth of the layer [M. On the whole, the heuristic 
Neel approach as well as similar models have succeeded to give qualitative explanations of some general features of the 
induced anisotropy. However, they do not account for many other observed effects and they fail to give a quantitative 
description of the magnetization structures in magnetic nanostructures. 

In this paper we propose a simple phenomenological theory which gives a consistent method to calculate the 
anisotropy and the equilibrium magnetization in thin films and multilayers. Let us consider a magnetic nanostructure 
embedded into nonmagnetic media. It may be a thin magnetic film confined by the interface with the substrate 
and the surface, or a magnetic layer between nonmagnetic spacers in a multilayer, or any other magnetic inclusion 
in a nonmagnetic matrix (such as nanowires, magnetic dots, or magnetic clusters). The magnetic energy of such a 
nanostructure can be written as an integral over its volume 
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Here m = M/Mo is the normalized value of the magnetization vector M (Mq = |M|), n is a unity vector perpendicular 
to the surface of the layer, A is the exchange stiffness constant, H is an external magnetic field, and is a 
demagnetizing field. The second term in (|^) is the energy density of the induced uniaxial anisotropy. To describe the 
distribution of K(r) within the nanostructure we introduce an interaction functional of Landau-Ginzburg type 
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The first term in (Q) represents the stiffness energy and f(K) is the energy of a homogeneous induced anisotropy. In 
many cases the function f(K) may be written as f(K) — aK 2 + 2bK. The stiffness parameter a > together with 
a > characterize the resistance of the system against the uniaxial anisotropy imposed by the boundaries. The linear 
term in f(K) reflects the tendency for symmetry breaking and the ensuing rise of homogeneous uniaxial anisotropy. 
Minimization of f(K) yields the value of this anisotropy Kb = —b/a . The scalar function K(r) may be treated 
as a physical field additional to the magnetization field M(r). Both fields are coupled by the anisotropy energy in 
(||). In most cases, this coupling energy should be negligible compared to the strong surface and volume interactions 
inducing the uniaxial anisotropy. Thus, the equilibrium distribution K(r) can be calculated independently from the 
magnetization field by minimization of the functional (||). The anisotropy on the confining surfaces is supposed to 
have values K(rjg — K Q fixing the boundary conditions for the variation problem. Then, the equilibrium distributions 
of the magnetization are determined by variation of the functional (|J) with this definite function K(v). Here, the 
boundary conditions Kq are kept constant. But, in general, Kq may vary within the surfaces, and all phenomcnological 
constants in (^|) may depend on characteristic sizes of the nanostructure. The values Kq (volume anisotropy density 
on the surfaces) should not be confused with the "surface anisotropy" Kg in (Q) which is an anisotropy energy per 
area. The relation between these quantities will be discussed below. 

As an application of the theory we calculate the induced anisotropy and the magnetization in a magnetic layer of 
thickness d sandwiched between two identical nonmagnetic spacers. The layer is supposed to be infinite in x and y 
directions and bounded by parallel planar surfaces z = ±d/2. On the surfaces the anisotropy is K(d/2) = K(—d/2) = 
Kq. The vector n is directed along the z-axis and Kq > (perpendicular, surface induced anisotropy). Within the 
magnetic plate K (r) varies along z-direction and is homogeneous in (x, y) plane. Variation of the functional (||) yields 
the following solution for K(z): 
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where A a = \J a/a is a characteristic length. To derive the equilibrium functions M(r) one has to minimize the 
functional for the magnetic energy (^) with K(z) given by (Q). We consider structures that are homogeneous in the 
plate plane (i.e. monodomain structures). Further, we ignore the variation of Mo observed in some multilayers [fl5|| . 
Under these assumptions the magnetization in the layer is described by the polar angle 6{z) between the vector M 
and z-axis. For such a one-dimensional problem the magnetostatic equations for have a rigorous solution, and the 
corresponding stray field energy is expressed by an anisotropy-like term Wd = — = 2ttMq cos 2 (8) [|16[ . This 

energy contribution is commonly named "shape anisotropy". Thus, the functional (|^) with K{z) given by (^|) may be 
expressed in the following reduced form 
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where h — H/ (2ttMq) is a reduced value of the magnetic field and ip is the angle between the vector H and z-axis, 
X ex = y/A/ (2ttMq) is the exchange length, and the volume density of the total anisotropy K(z) is 

K(z) - - 27 rM + cosh(d/(2Aa)) • (6) 

For simplicity we suppose that the homogeneous part of the induced anisotropy Kb (see (0)) is included into the 
coefficients Kq and 27rMg. 

The energy (|^) has the functional form of an uniaxial ferromagnet with anisotropy K(z) ^) varying along the layer 
thickness. The distribution of K(z) within the layer strongly depends on the layer thickness (Fig. 1). The length A a 
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characterizes the extension of the surface anisotropy into the layer depth and may be named "penetration depth of the 
induced anisotropy". When d X a the induced anisotropy is localized near surfaces and rapidly relaxes into the depth 
of the layer. For thinner layers, where A a and d have the same order of magnitude, the induced anisotropy attains 
considerable values within all the layer (Fig. 1). The possible magnetic phases and the regions of their stability are 
determined by optimization of the energy (p|) with free boundary conditions that read for this case (d6/dz) \ z= ±d/2 = 0- 
A complete analysis of the functional (g) will be done elsewhere. Here we only discuss some general properties of 
the solutions at zero field as representative for the problem. In this case the layer thickness d and the parameters 
q = Kq/ (27tMq) and A = A ex /A a span the phase space of the solutions. A typical (q, d) phase diagram for fixed 
values of A is shown in Fig. 2. It consists of two regions where homogeneous phases exist with either perpendicular 
(9 = 0) or parallel (9 = vr/2) magnetization. These are separated by an inhomogeneous (twisted) phase where the 
angle 9 gradually increases from a definite value 9i > on the surfaces to a largest value in the center 9\ < 9i < ir/2. 

The perpendicular phase is stable for small values of the thickness and when q > 1. The parallel phase exists 
only in rather thick layers and when the parameter q is smaller than a critical value q c (A). In the region q > q c the 
perpendicular phase transforms continuously into the twisted phase at the critical line di(q). For 1 < q < q c the 
region of the existence of the twisted phase is bound by two lines of second order transitions into the homogeneous 
phases di(q) and ^2(9)- Thus, according to our model the transition between the parallel and the perpendicular phases 
always occurs continuously via the transformation into an intermediary (twisted) phase. Fig. 3 shows the evolution 
of the magnetic structures during such a transition. 

By integrating K(z) in (^) over z one obtains for the average value of the uniaxial anisotropy density 

K 1 f d/2 w( W 9 n r2 , K tanh (d/(2A„)) 

Values of the anisotropy K e g have been measured in many layered systems Jl|,^] . Usually, experimental values of K e g 
are plotted as product K c g d versus d. In Fig. 4 the thickness dependence of the function $(d) = K e gd corresponding 
to (0) is shown for different values q. 

In the limit of large thickness (d^> X a ) the functions $(d) are negative and depend linearly on the thickness (with 
slope equal to the value of the shape anisotropy). In thiner layers (d < X a ) the functions #(<i) become nonlinear. For 
q < 1 they are always negative (parallel anisotropy) and reach zero as the thickness decreases to zero. For q > 1, the 
function <!>(d) changes sign at the critical thickness do determined by the equation 

and remains positive for < d < do manifesting the existence of the perpendicular anisotropy in this region. The 
function <P(d) reaches the maximum value for the thickness 
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and then monotonically decreases to zero as d — > (Fig. 4). Experimentally obtained functions <£(<i) are in close 
accordance with our theoretical results. E.g., nonmonotonic dependencies of (d) have been observed for Co/Au 
multilayers |T^| , Co/Ir multilayers or for Ni/Cu multilayers and Cu/Ni/Cu sandwiches p-|ll|. Taking experimental 
values for do and d max it is possible to calculate the parameters q and A a from eqs. (@) and (|j|). E.g., for Co/Au |l4| ] 
films this yields q — 4.6 and A a = 1.9A. For Ni/Cu multilayer systems the function <P(d) in the regions of perpendicular 
anisotropy has been obtained by three different experimental methods p|-pl|] . The results from these experiments show 
similar functional dependencies for <P(d). However, the measured values display considerable quantitative differences. 
Taking the experimental data from || and (l^] we derive the following values for the characteristic parameters q = 2.1 
and 2.2, A a = 26. 4A and 31. 9A. In contrast to the full curves from our theory in Fig. 4 experimental dependencies 
#(d) for these Ni-films become negative again for finite values of the thickness manifesting the reentrance of planar 
anisotropy in ultrathin films P-fij|- It should be noted in this respect that for the calculation of <P(d) (Fig. 4) the 
phenomenological parameters Ko are kept constant in this simple approach. According to experimental observations 
they may depend on the layer thickness, and they may be sensitive to details of the film preparation, e.g. via surface 
reconstructions, misfit strain or misfit dislocations etc. In particular, it was established experimentally that in the 
Ni-films under discussion surface anisotropy decreases with decrease of the layer thickness due to magnetoelastic 
interactions jL7|. For decreasing values of Kq the theoretical dependencies $(d) become negative in the limit of very 
thin films (dashed line in Fig. 4) . This thickness dependency of the phenomenological parameters Kg as well as other 
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phenomenological constants should be considered by introducing appropriate models in future work based on further 
experimental results and results from ab initio calculations. 

The distribution of the induced anisotropy and the magnetization in a magnetic layer can also be investigated 
within a discretized model which is useful in the limit of ultrathin layers or for comparison with ab initio calculations. 
In this case the equilibrium values of the induced anisotropy in the magnetic layer consisting of N magnetic planes 
are determined by minimization of the energy 
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where A is the distance between two adjacent planes. Equilibrium profiles of K n for symmetric boundary conditions 
A'i = Km = Kq are shown in Fig. 5. Using these profiles, the orientation of the magnetization vectors in the planes 6 n 
can be calculated by minimization of the discretized analogue of the energy (|J) . The characteristic functions for the 
discretized model <Pn = — 2itMqNA + XL=i K n shown in Fig. 6 display the same dependence on thickness d = NA 
as the continuous model (Fig. 4). 

Finally, we discuss the relation between our theory and the traditional approach. For large thickness (d ^> A a ), the 
function (||) may be asymptotically written as 



K{z) = X exp(- 
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The function K{z) has finite values only in the vicinity of the surfaces and equals zero in the main part of the layer. 
Integrating (O) with respect to z yields the anisotropy energy (per area) 
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The integral parameter K s is equal to the additional energy connected with surface induced anisotropy. Similar 
parameters characterize surface tension in liquids and domain wall energy in the theory of magnetic domains. Ex- 
panding K e ff in (Q) for large d/X a and substituting for K s (|lj) one obtains the Neel anisotropy ansatz ([!]). In the 
region of large thickness the surface anisotropy may be transformed into a surface term which should be included into 
the boundary conditions of the corresponding variation problem for the magnetization. For this case the calculation 
of the magnetization in magnetic layers firstly has been done by Thiaville and Fert Jl3| . It is clear that the surface 
anisotropy K s is meaningful only in the case of a strong confinement of the induced anisotropy to the surface, i.e. 
when d > A„. This relation determines the region of the validity for the Neel approach and other approaches that 
reduce surface effects to an effective surface energy contribution. Beyond this region the function K{z) of eq. (Eh has 
finite values within all the layer and the separation of the induced anisotropy into a surface and a volume contribution 
becomes meaningless. 

In conclusion, our theory considers inhomogeneous distributions of induced magnetic anisotropy. It is able to 
describe many of the experimental features found for thin ferromagnetic films. In particular, the nonlinear deviations 
from Neel's approach for surface anisotropy found in various experiments are consistently explained. 
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FIG. 1. Inhomogeneous distribution of the anisotropy K(z) across a magnetic film according to eq. (6). For the curves from 
the top to the bottom, the chosen value of A a decreases. 

FIG. 2. Typical phase diagram (schematic) for the magnetization structure of magnetic films described by eqs. (5) and (6) 
for fixed A = \ ex /\ a and in zero field. Depending on layer thickness d and on the strength of induced uniaxial anisotropy at 
the film boundaries q ~ Ko three phases may exist as symbolically depicted by the insets. A twisted phase is separated from 
the two homogeneous phases by second-order transition lines d\ and di. At the dashed line the energy for perpendicular and 
parallel phase would be equal. For q > q c only perpendicular and twisted phase exist. 

FIG. 3. Magnetization distribution across magnetic films for fixed q with increasing thickness d. Transition from perpen- 
dicular to twisted structure at di and from twisted to parallel at di- 

FIG. 4. Characteristic function 3>(d) = K c gd showing the dependence of the effective magnetic anisotropy eq. (7) on film 
thickness. Full lines give the typical behaviour when the anisotropy Ko at the boundaries is fixed. The dashed line shows 
reentrance of the parallel magnetization when Ko decreases with decreasing d in the case with q > 1. 



FIG. 5. Distribution of local anisotropy for the discretized model eq. (10) with various AT and fixed thickness. 

FIG. 6. Examples of characteristic functions for the effective anisotropy of the discretized model. Symmetric boundary 
conditions AT =0.05 (□), K =l (•) and \ 2 a = 10, b = 0.5, a = A = 1, M = 0.2. 
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